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Abstract 

Some results for the black-body radiation obtained in the context of the g-thermostatistics are 
analyzed on both thermodynamical and statistical-mechanical levels. Since the thermodynamic 
potentials can be expressed in terms of the Wright's special function an useful asymptotic expansion 
can be obtained. This allows the consideration of the problem away from the Boltzmann-Gibbs limit 
q = 1. The role of non-extensivity, q < 1, on the possible deviation from the Stefan-Boltzmann 
T 4 behavior is considered. The application of some approximation schemes widely used in the 
literature to analyze the cosmic radiation is discussed. 
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I. INTRODUCTION 



The nonextensive statistical mechanics (NSM) is based on the generalized entropy defined 

by: 



1 _ Y W v Q 



where the index i labels the possible microstates of the system under consideration, {p{\ is 
a set of normalized probabilities, the real parameter q characterizes the degree of nonexten- 
sivity and k is a positive constant. Notice that taking the limit q — * 1 leads to the popular 
Boltzmann-Gibbs statistics (BGS). For a recent review on the nonextensive thermostatistics 
and its current status see Refs. 0, 0| • 

The thermodynamics in the context of the (NSM) is investigated by generalizing the 
Gibbs canonical ensemble to the case q ^ 1. This is achieved by maximizing the entropy 
defined by Eq. (^) under the constrains: (i) normalization of the probabilities, and (ii) 
knowledge of the expectation value of the energy. The expectation values that lie in basis 
of the thermostatistical considerations are usually computed using two different approaches. 
The first one is the so the called 'unnormalized' approach proposed in Ref. Within this 
approach, for a given observable O with an eigenvalue Oj in the microstate % one has 

w 

<0> = J>?°«- (2) 

i=X 

This approach shows several difficulties in describing the thermodynamics (see e.g. (H). It is 
unable to preserve many of the thermodynamic properties. To overcome this inconveniences 
the 'normalized approach' has been advanced in Ref. where the expectation values are 
given by 

2-ii=\ Pi 

The normalized treatment seems to provide one with a natural bridge that connects the 
NSM to the thermodynamics ■ The normalized approach has been in turn improved by so 
called 'optimal Lagrange multipliers' (OLM) approach Nowadays it is believed that both 
approaches are the most appropriate choice for investigating the thermodynamics within 
the framework of the NSM regarding the nature of the Lagrange multiplier associated with 
the temperature. In these cases the nonextensivity is restricted just to the entropy, while 
the internal energy remains extensive as in the case of the BGS. The Lagrange multipliers 
preserve their traditional intensive character and can be identified with their thermodynamic 
counterparts. An important consequences of the success of both approaches is the unification 
of the Tsallis and Renyi variational formalism under a common context. This success is 
originating from the fact that the Renyi's entropy is extensive^. Let us mention that the 
extensivity of the internal energy is true as long as we consider a system with a large number 
of particles or the thermodynamic limit under the condition q < 1 Q . 

In the last few years many papers have been published on the application of NSM to 
the black-body radiation 0, H, M EH U, EH El; EH 15]. Our belief is that rigorous and 



(O) = (3) 



exact results have an instructive role in the field. The exact expression for the corresponding 
partition function has been obtained in Refs. [H EH 3- The studies presented in Refs. 
IE IH IH EH U) IH 13 employ the unnormalized approach, while those in Ref. [IH Hl| used 
the normalized one. 
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Irrespective of the used approach one can see that in the case q 7^ 1 the derivation of exact 
results is a very complicated task. The final expressions one has at hand are too cumbersome 
that obscures the underlying physics. In this situation an approximation that tries to make 
the exact results more simple and transparent is preferable. However, the real benefit from 
the exact treatment without well defined range of validity of the used approximation seems 
to be doubtful. Let us note that the more a ppr opriate ^approximations are limited to simply 



fsee Refs. 



computing (1 — q) corrections 
Boltzmann-Gibbs limit q — > 1 leads to great simplifications 



ll| and |l4| and references therein) since the 
In this situation the possible 
strong deviations from the usual Boltzmann-Gibbs case are of a significant interest. Quite 
recently the necessity to have a well estimated approximation in this field grows up, since 
some strong criticism concer ning the physical validity of Tsallis statistics takes place [la ll7T | . 
The objections of Refs. 3 12 1 i n their major parts are concentrated on applications of q- 
thermostatistics to the black-body radiation. Recently this issue has been a matter of a 
debate in the literature 3, 17. 18j|. 

The aim of the present study is to illustrate another possibility for simplification of the 
basic expressions in both approaches not related to the small value of (1 — q). It is based on 
the fact that in both cases the intricate sums that appear in the theory may be presented 
[l9l | in terms of the Wright function with well studied analytical properties [^j. This is 
justified since we are considering a tremendous system. We hope this possibility will shed 
some light on the existing debate 3, 17. 

This paper is organized as follows: In Section |H] we discuss the thermodynamic derivation 
of the popular Stefan-Boltzmann law in terms of the Tsallis statistics. In Section 11111 we 
introduce the mathematical background we need in our analysis. This is presented in Section 
II VI Section |V| is devoted to the discussion of our results. 



II. SOME THERMODYNAMIC RELATIONS 

First we shall introduce some basic notions. The radiation field in a large cavity can be 
considered to consist of a denumerably infinite set of electromagnetic oscillators correspond- 
ing to the various quantum states k in a <i-dimensional box. The oscillator frequencies, 
oji = cki, are related to the total energy E by E = YliUi^hjJi, where is the number of 
oscillator quanta with frequency oji and polarization e, c is the light speed, h is the Plank 
constant and ki = |k$|. The Boltzmann-Gibbs partition function Zi, for a large volume V, 
can be written as 

Z 1 (A d )=exp(A d ), (4) 

where 

A d= yMV±^( k -f) d V. (5) 



(47r) d / 2 r(d/2) V ^ c 

In (J5J), Trf = d— 1 is the number of linear-independent polarizations, /c^-Boltzmann constant, 
T-temperature, and T(x) and ((x) are Gamma and Zeta function, respectively. 

The main obstacles related to the applicability of the g-thermostatistics to the black- 
body radiation may be considered in the context of the famous Stefan-Boltzmann law. In 
the Boltzmann-Gibbs thermodynamics the Stefan-Boltzmann law follows from the equation 



U g (T, V) = kT 



and the relation 

P(T) = ^, (7) 

where p = p(T) is the pressure and u(T) = U (T, V)/V - the internal energy per unit volume. 
Here the dependence on the temperature alone is crucial. As a result Eq. © reduces to 
an ordinary differential equation for u(T) and its solution is u(T) = aT +1 , where a is a 
constant that cannot be obtained on the macroscopic level. 

Let us now consider the corresponding generalization of the Stefan-Boltzmann law to the 
NSM context. In the g-thermo dynamics the following expressions for the internal energy 
U q (T, V) and the pressure p g (T, V) hold 

2 d [Z q ]i-«-l 
dT (1-q) (8) 

and 

where Z q is the g-generalized partition function. Now, we shall give some thermodynamic 
relations for the black-body radiation using as input the definitions (J5J) and (jHJ). Because 
of the simple dimensional arguments it is evident that Z q = Z q (A d ). If we introduce the 
convenient notation ln g x = x<? 1 ^~ 1 the internal energy U q (T,V) may be expressed through 
Z q (A d ) and 

U q (T, V) = dkTA d -^- \n q Z q (A d ). (10) 
Correspondingly for the pressure p q (T, V) we get 

p q (T, V)V = kTA d -£- In, Z q (A d ). (11) 

From Eqs. (fTU|) and (fTT|) immediately follows the g-generalization of the relation (JJJ) 

P q (T,V)V=^^-. (12) 

The relation (|T2*|) between the pressure and internal energy is q- independent, as it should be. 
The violation of the relation between the pressure and the internal energy would compromise 
the theory since this can be established from pure electrodynamic reasoning. Eq. ()12j) was 
verified in Refs. and [l4| on the basis of the explicit expression of the partition function. 

A necessary and sufficient condition for the internal energy U q (T, V) to be proportional 
to the volume V and to obey Eq. (|1U|) is that Z q (A d ) must have the general form 

Z q (A d ) = e c q ^ A « (13) 

where e q — [1 + (1 — q)x] 1 ^ 1 ~ q ^ is the inverse function of \n q (x) and C\(q) is an unknown, 
regular at q = 1, function. Indeed the relation (|13j) cannot be exact. It can be obtained 
only as an approximation and C\{q) depend upon the used approximation scheme. For 
example, within the framework of the factorization approximation used in Ref. [lij . we 
have Ci(q) = [(4 — 3g)(3 — 2q)(2 — q)}^ 1 - Now, though g ^ 1 the Stefan-Boltzmann law 
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temperature behavior in its usual form is preserved. The constant o = a(q) must be in- 
dependent (see e.g. Refs. 0, B E3 Q ) • 

The same relation as Eq. © between p q (T, V) and U q (T, V) exists in the general case 
of the g-thermodynamics [21|. Here however instead of Eq. (J7J) the more general relation 
(112)) takes place and it is necessary to consider the following partial differential equation for 

This equation has a solution of the type 

U q (T, V) = aq (T d V)V c{q ' d)/d T 1+c{q ' d \ (15) 

where the constant C(q,d) (independent of T and V) and the unknown function u q {x) can 
be obtained only at the microscopic level. This is the generalization of the Stefan-Boltzmann 
law that can be obtained without using the explicit expression for the partition function. 

The result ()15|) means that in the considered case we loose the 'famous' T d+1 behavior 
of the internal energy as a function of the temperature. This is a strict consequence of the 
fact that U q (T, V) does not dependent linearly on the volume V. This is in agreement with 
the findings of Ref . 0] . However a consideration on a pure thermodynamic level does not 
exclude a g-dependence of the proportionality coefficient of the T 4 law. 

Without loss of generality let us consider a system in a cube with V = L d . If we introduce 
the mean thermal wavelength of the black-body photons I = l(T) = hc/kT, Eq. ()15)) may 
be transformed into the following scaling forms 

^i) = yg 9 (y), (16) 

where k is a dimensionless constant and g q (x) is a function, of which the explicit form 
depends on the way of writing the energy constraint (see Ref. 0]) i.e. its expression may 
be quite different as a function of the ratio L/l (see e.g. Eqs.(|22J) and (J2HJ)) depending on 
the used approach: unnormalized or normalized. 



III. THE WRIGHT FUNCTION 



In Section U we have advanced that for the investigation of the black - body radiation in 
the context of NSM different approaches has been used in the literature. Earlier, using the 
normalized approach the exact q counterpart of flU) is found to be [§] 



Zq ^ Ad) r (l -q) £ (1 -g) rfm m!r[(2-g)/(l - q) + dm}' ^ 



m=0 



Later, another expression for the partition function was obtained within the framework of 
the OLM and the normalized approaches. It is given 0] by the relation 

7(11 A\-T ( 2 -^A V A ? [1 + (1 ~ ^ kT ^} dm+1/{1 - q) 
q[q, d) ~ \l-q)^ {l-q) dm m\ T[(2 - q)/{\ - q) + dm] ' 1 ' 
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obtained under the cut-off-like condition 1 + (1 — q)(kT) 1 U q > 0, otherwise we have 
Z q (U q ,A d )=0. 

In spite of the fact that in the last case all the thermodynamic quantities, e.g. the internal 
energy, can be expressed in an exact fashion a complication arises. The corresponding 
expressions are self-referential [l3l . 14l | in the sense that the thermodynamic functions are 
not expressed in a closed form. This fact leads to mathematical difficulties that make the 
problem for the most part only numerically tractable. Notice that Eqs. (J 17)) and (|18|) are 
valid only for q < 1. 

If one tries to apply the above results to the experimental data of the cosmic back-ground 
radiation the condition Ad ^> 1 is always satisfied since hc/(kT) is of the order of 1/10 cm 
and V is of cosmological dimensions 0] . This physical fact will lead to a great simplification 
in the mathematical expressions given in Eqs. ()17)) and (fTHj) . Having this in mind, we take 
advantage of the fact that the series in the r.h.s of Eqs. (fT7j) and (fT8|) can be presented in 
terms of the entire function 



(f>(p,a;z) = 



m=0 



m\T(pm + ft) 



p > 0,a E C, 



(19) 



introduced in 1933 by E.M. Wright in the asymptotic theory of partitions. For analytical 
properties, some generalizations and applications of this function the interested reader may 
consult Ref. [2(j. We note here an useful mathematical result concerning the behavior of 
Wright function <p(p,a;z). If p > 0, for a large real z, we have the asymptotic expansion 

ET 



(1-2q) 

4>{p, ft; z) = (pz) ( 2 + 2 p) 



2tt 



exp^l+p- 1 )^) 



(i+p) 



M 

£ 

m=l 



a 4 p > a h o((pz)-m: 



(pz) ( 1+ "> 



(20) 



i.e. the asymptotic behavior of the Wright function is presented in terms of elementary 
functions. This result permits to obtain the different thermodynamic functions of the black- 
body radiation in a more simple form, namely, in some particular cases. The constants 
a m (p, a) can be exactly evaluated For our analysis below we need 



ai(p, a) 



and 



a 2 (p, a) 



P + l 



a 



ft 

2 {a 



p-l) + -(2 + p)(l + 2p) 



p) 2 L48 
7 



1152 



ft - p - l)[6ft 2 + ft(2 - Up) + p(6p - 7) - 2] 
(2 + p) 2 [103 + 4p(7 + p)] 



IV. THE STEFAN-BOLTZMANN LAW 

In order to define the unknown function and constants in the thermodynamic relations 
discussed in Section |TT] and to obtain the g-generalization of the Stefan-Boltzmann law one 
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must use the concrete expression for the partition functions: (fTTj) for Z q (Ad), obtained 
using the unnormalized approach, or ()18|) for Z q (U q , Ad), which is a result of the normalized 
approach. This motivates us to consider below both approaches separately. 



A. Unnormalized approach 



Within this approach the thermodynamic quantities are computed using the so-called 
unnormalized expectation values introduced in Eq. @. In this case the generalization of 
the Stefan-Boltzmann law in terms of the Wright function is given by the following expression 
for the internal energy 



U q (T,V) 



dkTA d 
(!-<?)< 



1 1-9 



d 



(1-qY 



d 2=2-- 

U > l-o> (1 



9' (1-9)° 



(21) 



Now, let us consider the physically interesting case d = 3. In the limit q — > 1 the 
asymptotic expansion (j2T?j) fails. Then if (1 — g) is fixed, for A 3 ^> 1, using Eq. (j2Hj) (up to 
the zeroth order in small values of z" 1 ) we get 



U g (T,V) 



kT 



( 87r )(l-„)/2 



r 



q 



1-9 



3 A, 



1 -(l-9)/8 



exp 



|[3A»(1 



?)] 



1/4 



(22) 



The last equation is in full consistency with the relation (|15|) if for the constant C(q,3) we 
take the value — 1(1 — q) and the function a q (T 3 V) is equal to the exponential function in 
the r.h.s. of Eq. (j22J) with the corresponding factor. 

In FIG. Q]we present the comparison of the behaviors of U q (T,V) from Eqs. (f2*Tj) and 
(|2*2*|) at large A3. It shows that for large A3, the expression (j21j) for the internal energy 
is a good approximation of the exact one given by (j22j) . To our knowledge such kind of 
approximations is presented for the first time for the black-body radiation problem, while 
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expansions around q = 1 investigating deviations from the BGS are known 

As a conclusion we find that our result (}2*2*|) shows that in the case of unnormalized 
approach the Stefan-Boltzmann law is not preserved. Remark that in this case the internal 
energy is not proportional to the volume. This is in concert with the discussion of Section 

El 



B. Normalized approach 

Within this framework the expectation values are computed using Eq. (JHJ). This approach 
has the advantage of reproducing the traditional thermodynamic relations. The internal 
energy U q = U q (T, V) is found to obey a nonlinear equation that can be expressed in terms 
of the Wright function. It has the form: 



4Ad[l + (1 " q^kT)- 1 ^]^ 1 ^{ d ^ q + d ' (1-9) 



U q = kT ai \ — qi — r-t-. (23) 

(l-9) d+1 6d,^-/ 



"> l- 9 ' {1-qY 
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FIG. 1: Behavior of g-dependence of U q (T,V)/3kT from Eqs. 1)21 Jl (solid line) and (1221) (dashed 
line) for (A) A3 = 500, (B) ^3 = 5000. Notice that we used the logarithmic scale along the vertical 



axis. 



For fixed (1 — q) and 



A d [l + (l-q)(kT)- l U q } d 

(i- q y 



>1, 



(24) 



using Eq. (|20j) (up to the first order in small values of z 1 ) after some algebra Eq. 
reads 



1 - q^kT)- 1 ^ = [1 + (1 - q^kT)' 1 ^} 

1 / dA d 
x 1 7 r 



[1 + (1 - q)(kT)~ l U q ] d 



(25) 
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The solution of this equation is surprisingly simple. The result is 



U q = aVT d+1 - -^—kT, (26) 

where 

r(d)C(d+l)2d(d-l) k d+1 
a ~ (4vr) rf / 2 r(d/2) Jhcj^ 

is the usual Stefan-Boltzmann constant. Since it is not possible to take the limit q — > 1 in 
Eq. f)26|) we cannot recover the Stefan-Boltzamnn any more. Inserting now Eq. ()26|) into 
Eq. (|24|) we obtain the condition 

A d {dA d ) d > 1, (27) 

restricting the range of the validity of the obtained solution for U q . In order to improve 
our result we calculated the next term of the internal energy f|26[). To this end we used the 
expansion of r.h.s of (|2*Hj) to the second term i.e. in small values of z~ 2 . The complicated 
ensuing equation could be solved using an iteration method, which leads to the additional 
term [1 — (1 — q)d\ [2d(l — q)}~ l kT in the solution (|2Tj) . Note here that in spite of taking into 
account the next order in our calculations we see that our results are not improving in the 
sense that we cannot take the limit q — > 1. 

Remark that the internal energy (J26|) is noextensive for relatively small volumes of the 
system. In the thermodynamic limit it becomes extensive, confirming the conclusions of Ref. 
p. The entropy, S, follows from the thermodynamic relation (dS/dV)x = (dp/dT) v . It 
has been demonstrated that in the case of the NSM the traditional thermodynamic relations 
remains valid only if one uses the Renyi entropy, S R , instead of its Tsallis counterpart 
22]. The first one has the remarkable property of being extensive. In the thermodynamic 



limit we have in the case under consideration 



S R = Va-^—T d . (28) 
d+ 1 v ' 

The Tsallis entropy can be deduced through the relation 

Sj=J^-(exp[(l-q)S R ]-l). (29) 
This is a sign that the Stefan-Boltzmann law remains valid in the NSM context as well. 



V. DISCUSSION 

Let us discuss the most important case d = 3. In our consideration a crucial point is the 
condition 

(30) 

that was used to truncate the asymptotic expansion ()2()j) for obtaining the results given by 
Eq. (El or Eq. flUJ). 

Our consideration shows that the application of the thermodynamical concepts of the 
NSM may lead to the T 4 Stefan-Boltzmann law. This takes place if the partition function 
has the form (fTB^) . This form would be a result of some approximations (see e.g. [lH) and 
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the T 4 behavior is a strict consequence of a linear on V dependence of U(T,V). What is 
important to note is that the inequality ()30|) prohibits the use of any formula of the type of 
Eq. (|13j) to the cosmic back- ground radiation. 

The criticism in Refs. [l6lll7l| against the use of the formula obtained under the condition 
that A 3 is considered as a small parameter and after that utilized to fit the data of the cosmic 
black-body radiation is supported by our consideration j23[. If one tries to apply the above 
results to the experimental data of the cosmic back-ground radiation the condition (J30|) is 
always satisfied as it was mentioned in Section UTTl since hc/(kT) is of the order of 1/10 cm 
and V is of cosmological dimensions. On the other hand, as it is mentioned fairly in the reply 
[l8| the criticism of Ref. concerns on equal ground both approaches: unnormalized and 
normalized which does not clarify the problem. In Ref. it is suggested that the last one 
is free of inconsistencies argued in Ref. 0] . In our understanding at this point the situation 
seems to be clear j24[. Here we, consider in more details both approaches unnormalized and 
normalized separately. 

In the case of the normalized approach our investigation results in the formula (|22j) . If 
the inequality (|30j) is fulfilled we loose the T 4 behavior and in addition the energy density, 
Eq. (1221). d epends on the volume, which remains unacceptable. This agrees with the results 
of Ref~Q. 

Free of such a defect would be a theory based on the normalized approach 0, 0, E] . 
In this case one can immediately see that the condition (|27p is the relaxed version of (|30|). 
This means that if one tries to apply the normalized approach to the analysis of the cosmic 
back-ground radiation the expression (|26|) has to be used. This result is consistent with 
the thermodynamic relations (|14|) and (|15|). Indeed the first term in Eq. (|26|) is the usual 
Stefan-Boltzman law. The question is: how to interpret the last one? The wisdom of the 
standard statistical mechanics is that such terms are to be omitted since they are of the 
order of 0(1/ V) and do not contribute in the thermodynamic limit. On the other hand, 
this term diverges with q — > 1~ and may be considered as a sign that in NSM the Boltzman- 
Gibbs limit q — > 1 and the thermodynamic limit do not commute with each other. This is 
in agreement with the results obtained in the framework of a classical gas j25| . In order to 
avoid any misunderstanding, let us emphasize that our treatment excludes any attempt to 
take the limit 1 — q and indeed U q > 0. 
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